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In this paper, we study a system of two-component Bose gas in an artificial magnetic field trapped 
by concentric harmonic and annular potentials, respectively. The system is realized by gases with 
two-internal states like the hyperhne states of 87 Rb. We are interested in effects of a Rabi oscillation 
between them. Two-component Bose Hubbard model is introduced to describe the system, and 
Gross-Pitaevskii equations are used to study the system. We first study the Bose gas system in the 
annular trap by varying the width of the annulus and strength of the magnetic field, in particular, 
we focus on the phase slip and superflow. Then we consider the coupled Bose gas system in a 
magnetic field. In a strong magnetic field, vortices form a Abrikosov triangular lattice in both 
Bose-Einstein condensates (BECs), and locations of vortices in the BECs correlate with each other 
by the Rabi coupling. However, as the strength of the Rabi coupling is increased, vortices start 
to vibrate around their equilibrium locations. As the strength is increased further, vortices in the 
harmonic trap start to move along the boundaries of the annulus. Finally for a large Rabi coupling, 
the BECs are destroyed. Based on our findings about the BEC in the annular trap, we discuss the 
origin of above mentioned phenomena. 

PACS numbers: 67.85.Hj, 74.25.Uv, 03.75.Nt 


I. INTRODUCTION 

In the past decade, ultracold atomic systems are one of 
the most actively studied fields in physics. It is expected 
that these systems can be a quantum simulator for vari¬ 
ous models in quantum physicsi. Another aim to study 
ultracold atoms is to search for new quantum phenomena 
whose existence is theoretically predicted. One of these 
example is the supersolid that has both a diagonal solid 
order and an off-diagonal superfluid order—. 

In this paper, we shall consider a system of two- 
component Bose gas in an artificial magnetic fields. Each 
component of the boson corresponds to an internal state 
of single boson. We consider a two-dimensional (2D) sys¬ 
tem, and one component of the boson, which we call A- 
atom, is trapped by a harmonic potential, and the other 
boson, called B- atom, by an annular trap. The atoms are 
coherently transferred by a Rabi coupling between two in¬ 
ternal states. It should be remarked that recently closely 
related atomic system was realized in experiments^. As 
we show in the present study, interplay of the harmonic 
and annular traps generates intriguing phenomena. That 
is, the system exhibits interesting behavior depending 
on the strength of the magnetic field, inter and intra¬ 
repulsions, and the Rabi coupling. The present study 
was originally motivated by the findings of co-existing 
two kinds of superconductivity in the Sn/Si core-shell 
clusters^, i.e., we expected that a phenomenon similar to 
the double superconductivity of the Sn/Si core-shell clus¬ 
ter can be realized in a ultracold atomic system. However 
in due course, we have realized that the ultracold atomic 
system has strong versatility and controllability. Much 
more interesting phenomena are to be observed in that 
system than we first expected. 

This paper is organized as follows. In Sec.II, we explain 


the system of the two-component Bose gas trapped by 
harmonic and annular potentials, and introduce Gross- 
Pitaevskii equations (GPEs) to study the system. The 
bosons are coupled with the artificial magnetic field that 
is generated, e.g., by rotating the system. The Rabi cou¬ 
pling that induces coherent transference between the A 
and 5-atoms is also introduced. Rather detailed descrip¬ 
tion of the numerical methods is also given. In Sec.Ill, we 
consider behavior of a single-component BEC in the har¬ 
monic and annular potentials separately. In particular, 
we study the ring-shaped BEC in the annular trap in de¬ 
tail. Behavior of the BEC is strongly depends on width of 
the annulus and strength of the artificial magnetic field. 
In the case of the quasi-ID case, a phenomenon similar to 
the Aharonov-Bohm (AB) effect takes place as recently 
observed by the experiments^. On the other hand for 
a relatively wide ring-shaped BEC, superflows are gen¬ 
erated in regions near the boundaries. The superflows 
are similar to the surface supercurrents in superconduc¬ 
tors in a strong magnetic field. In Sec.IV, we study the 
two-component Bose gas in the harmonic and annular 
traps with the repulsions and the Rabi coupling. We 
fix the intra and inter-repulsions, and vary the strength 
of the Rabi coupling to observe how the BECs change. 
For a finite but small Rabi coupling, vortices start to 
vibrate around their equilibrium locations. For a mod¬ 
erate Rabi coupling, vortices of A-condensate flow along 
the boundaries of the 5-atom ring-shaped condensate. 
As the Rabi coupling is increased furthermore, the BECs 
are destroyed. Section V is devoted for conclusion. 
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II. BECs IN HARMONIC AND ANNULAR 
TRAPS 

A. Potential traps and GPEs 

In this section, we shall derive GPEs describing two- 
component BECs trapped by concentric harmonic and 
annular potentials, respectively. The two-component 
bosons mean two-internal states of a single boson like 
the hyperfine states of 87 Rb and they can be trapped by 
different potentials by means of the difference in mag¬ 
netic/electric properties of thern^. We call these bosons 
A-atom and B-atom, respectively, and denote the state- 
selective potential by Va(b)( x )> where x = (aq,^) is the 
two-dimensional ( 2 D) coordinate. See FigHJ 



state-selective potential 



FIG. 1. (Color online) Harmonic and annular traps for two- 
component Bose gases. Potentials Va{b){ x ) are shown in the 
upper panel. The lower panel shows the superposed state- 
selective potential. 


When both the Rabi coupling g and inter-species repul¬ 
sion U are vanishing, the GPEs HUD are two independent 
equations, and the dynamics of the solutions are governed 
simply by the potentials Va(x), Vb(x) and the strength 
of the intra-repulsion R. As we show in the following 
sections, the inter-coupling V and the Rabi coupling g 
strongly influence the structure of the groundstate and 
low-energy excitations of the coupled system. In partic¬ 
ular, interplay of the Rabi coupling and the potentials 
Va(b) plays a very important role. 

We briefly explain the numerical methods to solve the 
GPEs (| 2 .ip . As the unit of the length, we used the spa¬ 
tial cutoff a = 0.5/im. Every length in this paper is mea¬ 
sured in this unit. For the energy unit, we used Na mass, 
mjv = 4.0 x 10 - 27 Kg, and then m ^ a2 = 80nK. Every 
energy is measured in this unit. Typical time slice used 
for solving the GPEs (j 2 .ip is r = 10 - 6 sec. The dissipa¬ 
tive damping factor 7 = 0 . 01 , as explained above, and we 
used the Crank-Nicolson method to evaluate ^a{b){PAt) 
from Vbi( b) (t) • The total numbers of particles were as¬ 
sumed 10 6 ^ 10 5 as in the typical experiments. In the 
initial state of the time evolution, the ratio of the to¬ 
tal number of atoms in the harmonic and annular traps, 
Ah and A a , was set Ah/A a = 1 . After each time evolu¬ 
tion, we normalized the wave functions 'ipA(B)it) to keep 
the total number of atoms constant. Typical number of 
the time steps is ^ 10 6 , and therefore typical time scale 
of numerical simulations is ^ lsec, which is comparable 
with the experiments. 

III. NUMERICAL RESULTS OF GPEs FOR 

SINGLE-COMPONENT BEC IN MAGNETIC 
FIELD 


GPEs of condensate-field ^a{x) and ips( x ) for the cou¬ 
pled A and B-atoms system are given by 

(* - = ( - & £i=i, 2 ($ - iA i ) 2 - M 

+ R^AtpA + UTpsi’B + V A jtpA + ffV’B, 

(i - = (-£ £*= 1 , 2 ^ - - » 

+ R'lpB'lpB + U^Ai’A + Vb)^B + g$A, (2-1) 

where R(U) is the intra(inter)-species repulsion and g is 
the Rabi coupling^. A{ is a vector potential for a uni¬ 
form artificial magnetic field, which is created by rotat¬ 
ing the system or using lasers. The parameter /jl is the 
chemical potential, which controls the particle number, 
and 7 is a dissipative parameter. We put 7 = 0.01 for 
most of the practical calculation. To solve the GPEs 
m, we employ the Crank-Nicolson method, and use 
the symmetric gauge for the vector potential A$, i.e., 
(Ai(x), A 2 (x)) = (tt/x 2 , — Trfxi). Therefore the mag¬ 
netic flux per unit area is 2nf. For a rotating system, 
/ oc mO, where U is an angular velocity of the rotation. 


A. Single-component BEC in artificial magnetic 
field 

In this section, we show the numerical results for the 
GPEs (|2.ip of the BECs in both the harmonic and an- 
nulus traps. We set ^ — 1. For the case of U = g = 0, 
which corresponds to the two independent BECs trapped 
by each potential, a stable vortex lattice is expected to 
form in each BEC for a finite intra-repulsion R > 0 and a 
strong magnetic field. We consider the case with the ex¬ 
ternal magnetic field / = 0.02, which is relatively strong, 
and R = 1. In FigJSJ snapshots are shown for the BEC 
confined in the harmonic potential. FigEJalso shows the 
time evolution of the BEC confined in the annular trap 
of the inner radius £ = 20 and the width h = 20, which 
is larger than the healing length of the BEC. To mea¬ 
sure vorticity v(x) of the BEC with an order parameter 
'ip = yfpe l6 , we use the following definition, 

v (x ) = 7 -—rot 6{x). (3.1) 

27r 

It is obvious that an Abrikosov lattice forms even for 
the BEC in the annular trap as verified by both the den- 
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FIG. 2. (Color online) The case of the vanishing Rabi coupling 
and inter-repulsion. Snapshots of density of the BECs in the 
harmonic and annular traps, and the locations of vortices v(x) 
in the annular BEC. Abrikosov vortex lattice forms in each 
BEC. Parameters are h — 20, — — 1 and / — 0.02. 


sity and vorticity snapshots in Fig]2l We have verified 
that as the strength of the magnetic field is increased, 
the density of vortex increases in both the harmonic and 
annular traps. The snapshots in Figj 2 ]clearly exhibit that 
BEC first becomes unstable by the sudden application of 
the external magnetic field and then subsequently vor¬ 


tices appear near the boundaries of the BEC, enter the 
bulk of the BEC and finally form an Abrikosov lattice^. 
Formation of the Abrikosov lattice of vortex in the annu¬ 
lar BEC implies the existence of stable surface superflow 
as we verified in the following section. 



FIG. 3. (Color online) Snapshots of BEC density in the an¬ 
nular trap with h — 10. Application of sudden magnetic field 
generates vortices but they move to the outside of the BEC. 
As a result, Abrikosov lattice does not form. Parameters are 
^ = 1 and / = 0.02. 

m J 


It is interesting to see how vortices behave in a nar¬ 
rower BEC confined in a quasi-one-dimensional (qlD) 
annulus. We studied the BEC in the annular trap with 
h = 10 , and show the numerical results in Figj 3 j Vor¬ 
tices are first generated by the sudden application of the 
artificial magnetic field with / = 0.02, but they move to 
the outside of the BEC, and no Abrikosov lattice forms. 
We also studied the BEC in the annular trap with h = 5 
and found that no vortices are generated by the exter¬ 
nal magnetic field. In the case with the magnetic field 
/ = 0.02 and h < 10, the boson system can be regarded 
as a ID system, and therefore it is expected that the su¬ 
perflow behaves differently from that in the 2 D case with 
larger h. In the following section, we shall verify this 
expectation. 


B. BEC in annular trap, weak link and phase slip 



FIG. 4. (Color online) Annular trap with a weak link. Poten¬ 
tial for the weak link is explicitly given by Eq. iET2ti . 
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The results in the previous section were obtained for 
the magnetic field / = 0 . 02 , for which the corresponding 
magnetic flux in the central hole region of the annulus, 
4>i n , is about Sir and very large. As a result, the vortex 
lattice forms in the bulk of the ring-shaped condensate 
with the width h > 20. On the other hand in most of the 
practical experiments on the rotating BEC in an annular 
trap, the total magnetic flux 4>i n is about tt ~ 47 r, and 
a phenomenon similar to the AB effect was observed. In 
some experiments, a weak link (WL) was produced in 
condensates by a focused blue-detuned laser field and a 
rotating BEC was produced by stirring the condensate 
by the WL. Recently, a specific method was invented 
and used to measure the phase of the ring-shaped BEC 
directly^. 

In order to describe the above mentioned experimental 
setup, we add a WL potential V w \ in addition to Vb(x). 
In the rotating frame, the WL potential is given as 

V w \ = w6(x i + 5)6(5 - xi)0(x2 — t)0(£ + h — X 2 ), (3.2) 

where w is the height of the potential barrier, 5 is its 
width, and we typically take <5 = 5. See Fig0J 

We expect that as the width of the annulus h is varied, 
behavior of the annular condensate changes gradually. In 
particular, in the sufficiently small h compared to the 
magnetic length determined by / and also the healing 
length of the BEC, the system is effectively described 
by a ID model. Corresponding to the experiments, we 
consider the case such that 4>i n ~ 2tt (i.e., / = 0.005). In 
this case, it is expected that no vortices exist in the bulk 
of the condensate even for a wide annulus like h ~ 20 . 

We first study the ring-shaped condensate with h = 5 
and — = 10. We verified that for w = 0 and f = 0 (i.e., a 
homogeneous condensate in the vanishing magnetic field) 
the condensate has a uniform phase, as is expected. 

For the numerical study, we first prepared the homoge¬ 
neous condensate for h = 5, re = 0 and / = 0, and then 
we suddenly applied the WL of w = 10(100) and the 
magnetic field / = 0.005. Solution of the GP equation 
shows that the condensate was first disturbed by the ap¬ 
plied magnetic field and WL, but became stable by the 
effect of the dissipation. The density and phase of the 
final states are shown in Fig{5l 

The results in Figj5j in particular the phase of the con¬ 
densate, should be compared with the experimental ob¬ 
servation. As reported in Ref . 7 for the rotating BEC, the 
phase of the condensate exhibits a sudden increase in the 
WL. As the superfluid (SF) velocity v s in the rotating 
frame is given by v s = ^W(x), the above behavior of 
6(x) means that SF in the WL flows in the opposite di¬ 
rection against to the bulk superflow. This result is in 
agreement with the experimental observation in Refs^, 
and theoretical studies in Ref a 1 Q i n . However, the calcu¬ 
lation in Fig J5] also shows the existence of smooth increase 
of the phase in rather wide region of the opposite side of 
the WL. This state in Figj5] is the lowest-energy state, 
and the result indicates that a phase increase in the WL 
induces another phase increase in its opposite side of the 



FIG. 5. (Color online) Density p and phase 6 of the ring- 
shaped condensate with WL. Height of the WL, w = 10 (up¬ 
per panel) and w — 100 (lower panel), / — 0.005 and — 10 . 
The phase of BEC increases in the WL and also the opposite 
side of the WL. 



FIG. 6. (Color online) From the states in Fig El we decreased 
/ and w. Density and phase of the ring-shaped condensate 
without the magnetic field and WL. It is observed that the 
state with the unit winding is realized. 


ring-shaped condensate. The supercurrent j observed in 
the static frame , j = ^(V — A)Q, on the other hand, ex¬ 
hibits no anomalous behavior as we show in the following 
section. 

Finally from the state shown in Figj5j we gradually 
decreased the magnetic field and WL to / = w = 0, 
and then observed the density and phase of the conden¬ 
sate. The results are shown in Figj6l This is obvious 
that homogeneous superflow is generated with the unit 
circulation^. We also studied the case of wider annular 
trap. The results are shown in FigJ3 For large w = 100, 
the phase of the condensate is almost constant except in 
the WL region. This result means that the BEC moves 
with the WL, as is expected. 

For a small hopping case like ^ = 0.1, the obtained 
results are shown in FiglHJ It is obvious that the deviation 
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FIG. 7. (Color online) Density and phase of the ring-shaped 
condensate with WL. Height of the WL w = 0,10 and 100, 


h — 20, / — 0.005 and — — 10. 

’ J m 



FIG. 8. (Color online) Density and phase of the ring-shaped 
condensate with the WL. Height of the WL w — 10 and 100, 


f = 0.005 and ^ = 0.1. 

m 


h 2 

of the phase is small compared to the case with = 10. 
For small hopping amplitude, the condensate moves with 
almost the same velocity with that of the WL, i.e., the 
SF velocity in the rotating frame is very small. 


C. Supercurrents 

In this section, we shall study the supercurrent j (x) in 
the ring-shaped condensate. In the previous sections, we 
showed that there are two distinct states of the BEC in 
an external magnetic field, i.e., in one case vortices exist 
in bulk of the condensate because of a relatively large 
width h and also a strong magnetic field. The other case 
corresponds to the thiner condensate and/or condensate 
in a weak magnetic field in which no bulk vortices exist. 
Different behavior of the supercurrent is expected for the 
above two cases. 

We first consider the case of the ring-shaped conden¬ 
sate with / = 0.02, h = 5 and the WL w = 10. The 
result is exhibited in Figj9j which shows the supercur¬ 
rent j in the azimuthal direction. It is obvious that the 
supercurrent flows clockwise. Near the WL, the current 



-50 0 50 X 1 


j(x) 



FIG. 9. (Color online) Supercurrent of the ring-shaped con- 

h 2 

densation in the static frame. Parameters are — = 10, 

m 7 


/ a 0.02, h = 5 and w m 10. 



J [x. 



FIG. 10. (Color online) Supercurrent of the ring-shaped con- 

j- 2 

densation in the static frame. Parameters are — = 10, 

m ’ 


/ = 0.02, h = 20 and w = 10. 


is vanishingly small because of the low density of the con¬ 
densate. As the supercurrent is conserved in the state of 
a static density, there exists supercurrent in the radial 
direction near the WL as the previous study showedii. 

On the other hand for the BEC in which vortices exist 
in the bulk, the supercurrent flows as shown in FigfTOl 
From the result, it is obvious that the there exist sur¬ 
face supercurrents near the inner and outer boundaries 
and they are stabilized by the vortex lattice in the bulk. 
This behavior resembles the ordinary second-kind super¬ 
conductor in a strong external magnetic field. As in the 
quasi-one-dimensional case in Figj9j the supercurrent has 
to have a radial component near the WL. This surface 
supercurrent plays an important role when two BECs in 
the harmonic and annular traps interact with each other 
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through the Rabi coupling, as we show in Sec.IV. 


D. Inter-species repulsion vs. Rabi coupling 

By the practical calculation, we verified that the Rabi 
coupling tends to make the A and 5-atom condensates 
overlap with each other. On the other hand, the inter¬ 
species repulsion makes two condensates immiscible. In 
this section, we study how the overlap of two conden¬ 
sates varies as the ratio of the Rabi coupling and the 
inter-species repulsion, changes. To this end, we trap 
both the A and 5-atoms in the same harmonic type po¬ 
tential. The depth of the potential is controlled to realize 
a substantial overlap of the two condensates, i.e., we em¬ 
ploy a deep harmonic trap. In the practical calculation, 
we put U =1.1, 5 = 1.0, and vary the Rabi coupling g 
to observe how the condensates change. The results of 
the calculation are shown in FigJTTJ 



(c) g=0.03 


(d) g=0.05 


FIG. 11. (Color online) Density overlap of A and 5-atoms’ 
condensates. Each BEC contains three vortices. As the Rabi 
coupling is increased, the larger overlap appears. The inter¬ 
species repulsion U = 1.1. 


In FigQTj as a result of the strong magnetic field, there 
exist vortices in each condensate. For the case of weak 
Rabi coupling as g = 0.01 and g = 0.02, vortex cores 
of the two condensates are separated. As g is increased, 
the locations of cores are getting closer, and for g > 0.05 
they coincide. 

The above observation is useful to understand proper¬ 
ties of the condensates trapped in two different but finite- 
overlap potentials, which are investigated in the following 
section. 


IV. BECs IN HARMONIC AND ANNULAR 
TRAPS: RABI COUPLING AND INSTABILITY 
OF BECs 


In the present section, we shall study dynamics of two 
BECs trapped by the harmonic and annular potentials, 
respectively. We call the atom trapped the harmonic (an¬ 
nular) potential A-atom (5-atom) as before, and these 
atoms are coherently transferred by the Rabi coupling. 
In the practical calculation, we fixed the inter and intra¬ 
species repulsions U and 5, and then we varied the Rabi 
coupling g as a controllable parameter. At t = 0, we sud¬ 
denly add the Rabi coupling and then we simulate how 
the condensates evolve. We show the results for the case 
with / = 0.02, U = 0 and R = 3.0, although similar 
results were obtained in the case of U >0. 





t= 100 msec 









[ y .1 





t 

= 200 msec 


t 

= 300 msec 


FIG. 12. (Color online) Density profiles of condensates in 
the harmonic and annular traps as a function of time. Rabi 
coupling g — 0.01 is added at t = 0. Both BECs reach stable 
states in which stable vortex solids form and vortex locations 
are correlated with each other. / = 0.02. 


Ja(x) 

Jb{x) 




FIG. 13. (Color online) Time evolution of the supercurrents 
after a sudden turning on the Rabi coupling g = 0.01. The 
existence of the Rabi coupling does not substantially influence 
behavior of the supercurrents. 









































































7 


We first show the results for the case of a weak Rabi 
coupling, i.e., g = 0.01. See FigJT2] for time evolution 
of the condensates. After a sudden turning on the Rabi 
coupling, both the BECs form stable states in which vor¬ 
tex solid is stable. We also show the behavior of the 
supercurrents in FigUHl It is obvious that the existence 
of the Rabi coupling does not substantially influence be¬ 
havior of the super currents. The ring-shaped condensate 
exhibits surface supercurrents near the inner and outer 
boundaries, whereas the condensate in the harmonic trap 
has no specific pattern of the supercurrent although the 
supercurrent has a strong correlation with the vortex lo¬ 
cations in the bulk of the condensate. 




FIG. 14. (Color online) Beahvior of vortices after a sudden 
turning on the Rabbi coupling g = 0.05. Vortices vibrate. 


1a(x) 

3b(x) 



-50 0 50 -50 0 50 

X\ X2 

X 2 = 0 t=0 msec x\ = 0 



FIG. 15. (Color online) Supercurrents at t = 0 msec and 
300 msec. At t — 300 msec, it is observed that the A and 
B supercurrents synchronize with each other. Rabi coupling 
g = 0.05. 


Next let us consider the case g = 0.05. By the prac¬ 
tical calculation, we observed time evolution of the con¬ 
densates, in particular, we were interested in behavior of 


FIG. 16. (Color online) Vortices of A-atom condensate move 
along the boundaries of the ring-shaped B-atom conden¬ 
sate. In the outer (inner) boundary, vortices move clockwise 
(counter-clockwise). Rabi coupling g = 0.1. 



FIG. 17. (Color online) Supercurrent generated in A-atom 
condensate by the Rabi coupling. Direct movement of A- 
atom vortices induces the clockwise and counter-clockwise 
supercurrents that are correlated to the supercurrents in the 
B-atom condensate. Rabi coupling g — 0.1. 


vortices. See Fig(14]for the numerical results. We found 
that vortices vibrate around their equilibrium locations. 
Precise investigation shows that vibration of A-atom vor¬ 
tices is larger than that of B-atoms. By the Rabi cou¬ 
pling, transference of A and B-atoms takes place, and as 
a result, an effective interaction between A-vortex and 
B-vortex appears. The repulsion between intra-species 
vortex and the effective inter-species vortex attraction 
by the Rabi coupling compete and they cause the vor¬ 
tex vibration. In Figll5l we show the supercurrents in 
the condensates at t = 0 and t = 300 msec. At t = 300 
msec, it is observed that the A and B-atom supercurrents 
partly synchronize. 

For the case of g = 0.1, we found that vortices of A- 
atom move along the boundaries of the ring-shaped B- 
atom condensate. See Figjl6l In the outer (inner) bound¬ 
ary, A-vortices move clockwise (counter-clockwise). This 
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t = 0 msec t = 100 msec 



t = 200 msec t = 300 msec 




t = 200 msec t = 300 msec t = 200 msec t = 300 msec 



FIG. 18. (Color online) Densities of the condensates and su¬ 
percurrents in the annular region of both A and 5-atom con¬ 
densates. Rabi coupling g — 0.25. As a result of relatively 
large Rabbi coupling, A atoms concentrates in the region of 
the annular trap. 


FIG. 19. (Color online) Densities of the condensates and su¬ 
percurrents in the annular region of both A and 5-atom con¬ 
densates. Rabbi coupling g = 1.0. As a result of large Rabbi 
coupling, A atoms concentrates in the region of the annu¬ 
lar trap, and 5 atoms have a substantial density only near 
the boundaries of the annulus. No vortices are observed and 
there exit no stable supercurrent. This result indicates that 
coherent BECs become unstable by the large Rabbi coupling. 


is a result of the Rabi coupling and the surface supercur¬ 
rents in the 5-atom condensate. See also supercurrent 
in the A-atom condensate generated by the Rabi cou¬ 
pling in FigJlTl It can be seen that the supercurrents 
in the 5-atom condensate induce the ring-shaped super¬ 
currents in the A-atom condensate. Direct movement 
of A-atom vortices generates the clockwise and counter¬ 
clockwise superflows correlated to those in the 5-atom 
condensate. 

We furthermore increased the Rabi coupling g to g = 
0.25. Densities of the condensates and supercurrents are 
shown in FigHHJ As a result of relatively large Rabi cou¬ 
pling, A-atom tends to concentrate in the region of the 
annular trap, and the density of A-atom varies as a func¬ 
tion of time. On the other hand, the 5-atom condensate 
seems relatively stable. It is very interesting to notice 
that the A-atom condensate has a supercurrent in the 
ring azimuthal direction and it flows only in one direction, 
i.e., clockwise, whereas the 5-atom condensate has sur¬ 
face supercurrents flowing clockwise and anticlockwise, 
respectively. 

Finally, we show the results for the strong Rabi cou¬ 
pling g m 1. See Fig HU A-atoms concentrate in the 
region of the annular trap, and 5-atoms have a large 
density near the boundaries of the annulus. Vortices dis¬ 
appear and no stable supercurrents exist. This result 
means that coherent BECs become unstable by the large 
Rabi coupling. 


V. CONCLUSION 

In this paper, we studied the two-component Bose gas 
trapped by the harmonic and annular potentials. They 
interacts with each other by the coherent Rabi coupling 
and the inter-repulsions. The GPEs were used for study¬ 
ing low-energy states of the system. 

We first investigated single-component bosons in the 
harmonic and annular traps, respectively. In particular, 
the detailed study on the ring-shaped condensate in the 
annular trap showed that vortices exhibit interesting be¬ 
havior by applying an artificial magnetic field. We also 
studied the effect of the WL located in the annular con¬ 
densate, and obtained the phase of the BEC by solving 
GPE. Compared with the previous calculation in the ID 
case, we found the unexpected behavior of the phase of 
the condensate with the WL. We measured the super¬ 
currents and showed that they exhibit various different 
behaviors depending on the width of the annular conden¬ 
sate. 

Finally, we investigated behavior of the BECs in the 
coupled system with the Rabi coupling and repulsions 
between the A and 5-atoms. As the strength of the 
Rabi coupling is increased, the Bose gas system exhibits 
various “phase transition”. First A and 5-atom vortices 
start to correlate and they vibrate around their equilib- 
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rium locations of the vortex lattice. As the Rabi cou¬ 
pling is increased, vortices of the A-atom move along the 
boundary of the ring-shaped 5-condensate. This phe¬ 
nomenon occurs as a result of the superflow at the bound¬ 
aries of the 5-atom BEC and the Rabi coupling. This 
observation indicates an interesting possible method to 
control vortex movement in the condensate. We found 
that for a very large Rabi coupling BECs become unsta¬ 
ble as the vortex flow disturbs the A-atom BEC and then 
it induces instability of the 5-atom BEC again via the 
large Rabi coupling. 

The results in the present study have been obtained 
by solving the GPEs. It is desirable to see how the fluc¬ 


tuations, which are not taken into account in the GPEs, 
change the results. The path-integral Monte-Carlo (MC) 
simulations are useful for this study, although the time 
evolution of the system cannot be investigated by the MC 
simulations. The MC simulations on the present system 
are under study, and we hope that the results will be 
reported in the near future. 
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